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Abstract 

We consider polymers formed from a (quenched) random sequence of charged 
monomers of opposite signs. Such polymers, known as polyampholytes (PAs), 
are compact when completely neutral and expanded when highly charged. We 
examine the transition between the two regimes by Monte Carlo simulations, 
and by analogies to charged drops. We find that the overall excess charge, 
Q, is the main determinant of the size of the PA. A polymer composed of N 
charges of ±qo is compact for Q < Q c k, q ^/~N, and expanded otherwise. The 
transition is reminiscent of the Rayleigh shape instability of a charged drop. 
A uniform excess charge causes the breakup of a fluid drop, and stretches 
out a polymer to a necklace shape. The inhomogeneities in charge distort the 
shape away from an ordered necklace. 
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I. INTRODUCTION 



Polyampholytes (PAs) are long chain macromolecules with a random mixture of oppo- 
sitely charged groups fixed along their backbone IJ. Several elements conspire to make 
the behavior of such heteropolymers a problem of great interest. One (admittedly some- 
what remote) motivation is the similarity to the macromolecules of biological interest such 
as nucleic or amino acids. The specific sequence of monomers on such chains is essential 
to biological activity. For example, the sequence of amino acids determines the ultimate 
shape of a protein 0. Attempting to unravel the precise factors responsible for protein fold- 
ing, several statistical models have been proposed ||. These models sacrifice the specifity 
of particular proteins, by essentially focusing on generic properties of heteropolymers with 
competing interactions M. PAs can be regarded as a particular example of this class. From 
another perspective, properties of random systems with competing interactions has been on 
the forefront of statistical mechanics for the past decade The prototype of complexity 
in this class of problem, the spin glass, has much in common with random heteropolymers. 
The statistics of the ground state, and those of low lying excitations, is paramount to both 
systems. As examples of soft condensed matter, heteropolymers have the advantage of faster 
relaxation, compared to their "harder" counterparts. Indeed there is much encouragement 
from recent experimental studies of solutions and gels composed of PAs 

It may appear that the long range nature of the Coulomb interaction between charges is 
yet an additional complication of an already hard problem. Yet for a uniformly charged a 
polymer (a polyelectrolyte) , it is possible to find the exact scaling of the radius of gyration 
R g on the number of monomers iV |§. The proof relies on the non-renormalization of 
the dimensionless interaction parameter u = Q 2 /{ksTRg 1-2 ) in d embedding dimensions. 
Since Q oc N, it follows that R g ~ N u , with v = u h = 2/{d - 2) (for 4 < d < 6). 
Inspired by the simplicity of the homogeneous case, we suggested a similar argument for the 
randomly charged PA ||10|| . Consider a model PA composed of N charges ±go, randomly and 
independently chosen at each site. Although the mean net charge is zero, a typical PA has 
excess charge of order ±q ^/N. Independent of its sign, this leads to a repulsive self energy, 



on average described by the dimensionless parameter u = Nq 2 /(kBTRg ). Assuming that 
this parameter is not renormalized as in the uniformly charged case leads to a swelling 
exponent of v — l/(d — 2) (for 3 < d < 4), i.e. a polymer that is stretched in d = 3. 

However, in an electrolytic solution oppositely charged ions rearrange so as to screen 
the long range Coulomb interaction. The net effect is an attractive energy [[l l[, described 
by the Debye-Huckel (DH) theory. Higgs and Joanny [0, assumed that the monomers in 
a PA can similarly rearrange to compact configurations, thereby taking advantage of the 
DH attraction. A partial resolution of the contradiction between the two predictions is 
obtained by noting that DH theory requires the exact neutrality of the electrolyte, while 
the RG-inspired approach depends on the excess charge in a typical sequence. Monte Carlo 
simulations [0,0] indeed confirm that PAs with Q = compactify at low temperatures. 
By contrast sampling all random quenches with unrestricted Q produces a broad range of 
sizes, with an average consistent with R g oc N. 

Depending on the conditions, experiments observe both compact and expanded confor- 
mations 0. One set of experiments [JF] were performed on gels produced by crosslinking 
PAs ||. By changing the conditions of the solvent (pH, salt content, etc.), it is possible 
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to control both the excess charge Q on the PA chains, and the screening length. Due to 
the large screening length, the Coulomb interactions are significant. As a function of Q, 
the gel undergoes dramatic changes in volume; the neutral gel being the most compact. 
However, the volume of the gel does not change gradually with increasing charge. There 
is an interval of Q around the neutral point where the gel remains compact, suddenly in- 
creasing in volume by a order of magnitude beyond a threshold Q c . The threshold scales 
with the number of monomers N, within a screening length as Q c ~ qo\/N. Motivated by 
the experimental results, we undertook a systematic examination of the dependence of the 
size of a PA on its excess charge. Monte Carlo simulations show that PAs are compact for 
small Q, and expanded when Q exceeds a critical value of Q c ~ qoy/N, in complete analogy 
with the experiments. Some aspects of this transition can be understood by analogy to the 
behavior of a charged drop. The spherical drop is stable for charges smaller than a Rayleigh 
limit Qr oc y/V, where V is the volume of the drop. Charged beyond this point, the drop 
elongates to minimize the Coulomb repulsion. The elongated drop rapidly disintegrates into 
smaller droplets. In attempting to follow a similar scenario, a PA chain breaks into a neck- 
lace of globules connected by strings. The detailed shape of the necklace is determined by 
charge inhomogeneities. 

The main results of this work were summarized in an earlier publication [T5 |. Here we 



provide more detailed results and additional information. The paper is organized as follows: 
The competing arguments applied to PAs are discussed in some detail in Section [0], and their 
inconsistencies are emphasized. The results of Monte Carlo simulations on the dependence 
of the radius of gyration of the polymer on temperature and excess charge are presented in 
Section [II]. Various details of the Monte Carlo procedure are relegated to Appendix [A[ We 
argue that as Q is increased beyond a threshold Q c , the PA undergoes a sudden transition 
from a compact to a strongly elongated state. In Section |TV] we provide a qualitative picture 
of this transition by analogy to the shape instability of a charged drop. Some known results 
pertaining to such drops, as well as new calculations for spheroidal shapes are presented in 
Appendix |B[ Such analogies cannot be extended to the strongly distorted limit where, as 
argued in Section a uniformly charged polymer deforms into a necklace of compact beads. 
Such a shape is the best compromise for the PA in trying to mimic the ground state of a 
charged drop, which as described in the Appendix |C], is obtained by splitting into several 
droplets. In Section [VT] we point out the importance of quenched randomness in the PA. An 
ordered necklace is not stable to charge inhomogeneities, and the beads must rearrange in 
complicated shapes dependent on the details of randomness. 



II. POLYAMPHOLYTE PHENOMENOLOGY 

In this work we consider PAs immersed in a good solvent in which the concentration 
of counterions is small, and hence the electrostatic interactions are treated as unscreened. 
Experimentally, the details of a charge sequence are determined by its fabrication process. 
Markovian sequences are constructed by adding one monomer at a time, the probability of 
choosing a particular monomer (e.g. positively or negatively charged) depends only on the 
last monomer fTE| . Such a construction leads to correlations in the charges g,, which decay 
exponentially as, 



3 



,2 



A |i_il , -1 < A < 1. (1) 



= 9o- 

The extreme limits of A = —1 and A = 1 correspond to non-random sequences that are 



alternating [17|, or fully charged, and will not be considered here. The behavior for A is 
asymptotically similar to the case where the charges are uncorrelated (A = 0), and we shall 
focus mostly on such chains. 

We shall, however, allow for the possibility that the total charge of the chain is con- 
strained to a particular value of Q. Experimentally, this net charge can be controlled by 
changing the pH, and other properties of the solvent. At several points we shall also contrast 
the behavior of quenched and annealed PAs; the latter is defined a predetermined collection 
of positive and negative charges which can freely move along the polymer chain. Note that in 
this definition the total charge is fixed. If this constraint is also removed the polymer lowers 
its energy by getting rid of any excess charge. Note that such annealed PAs are not very 
realistic, and should not be confused with the case considered by Raphael and Joanny 



where the excess charge can adjust its value depending on the temperature, concentration 
of PAs and concentration of the counterions. 

A short ^-monomer segment of a PA with uncorrelated random monomers has typical 
charge of qo\fl. If we assume that the segment is a self-avoiding walk, its radius of gyra- 
tion is approximately a£ u , where a is a microscopic length (e.g. the monomer diameter or 
nearest neighbor separation along the chain). As the typical electrostatic energy of such a 
subchain is q$£ 1 ~ u /a, interactions become important for T « q\£}~ v ja. (We shall hence- 
forth measure temperature in energy units, i.e. set ks = 1.) Alternatively, we can define 
£t = (Ta/qQ) 1 '( 1 ~ v \ and divide the entire chain into segments of It monomers. The in- 
teractions within each segment are small compared with T, while interactions between the 
segments are strong. Such segments form the basic "ions" in an analogy |II§ between PAs 



and usual electrolytes |! !| employed by Higgs and Joanny ||12|| . The spatial extent of each 
segment is 

a T = at T = I a ( f ) Vl{1 ' V) > for T > til a ' (2) 
[a , for T < q^/a. 

For the generalized Markovian chains, the net charge of a sufficiently long segment (such 



that A <C 1) grows as qoJN(l + A)/(l — A). The previous argument thus remains valid 



after replacing q by goy (1 + A)/(l — A). This simple change is sufficient to relate most 
high temperature properties of the short-range correlated sequences to the uncorrelated 
ones (e.g. in a high temperature series expansion for the radius). 

In the spirit of an RG analysis, we can attempt to increase the short distance cutoff along 
the sequence by a factor I. Perturbations around the high temperature phase depend on the 
dimensionless interaction parameter u = q}/r{r <C 1, where qe is a typical charge on length- 
scale £, while re is the spatial extent of the shortest segments. Upon rescaling the cutoff by 
a factor b, qt increases by \fb, while re scales by b v , as in a self-avoiding walk. Thus, the 
renormalized interaction parameter, u(b) = b x ~ v u grows under rescaling and reaches unity 
for £ equal to £t defined above. At this scale re becomes of the same order as ax- Beyond 
this point the interactions are relevant, strongly modifying the behavior of the chain. 

Different approaches to the problem are in agreement in the weak coupling regime of 
ii< 1. The strong coupling regime is not easily tractable, and different assumptions lead 
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to different conclusions. Higgs and Joanny [12| construct the free energy of a long PA 



by first assuming that it has a uniform density and then proving this assumption self- 
consistently. (See also Refs. In this approach the iV-monomer PA is divided into 



blobs |2(| of It monomers each, forming a liquid of uniform density as depicted qualitatively 
in Fig. |l]a. The blobs are non-interpenetrating and arranged so that the neighborhood of 
each blob is predominantly occupied by blobs of opposite charge. This arrangement roughly 
resembles the structure of a salt crystal. The excess charge of each blob is effectively screened 
and the dense configurations take energetic advantage of the large number of neighbors of 
opposite sign. The energy gain per blob is approximately the nearest neighbor interaction, 
i.e. e c (T) w q\£ T jaT- We note, however, that even in the ideal NaCl crystal, the condensation 
energy per atom (0.874gp/a) is a small fraction of the interaction energy between the nearest 
neighbors (3gg/a), and is thus strongly influenced by further neighbors. The validity of the 
picture depicted in Fig. |l]a rests on the assumption that it is possible to fold a randomly 
charged object in a way that not only provides the correct neighborhood to each charge, but 
also keeps more extended neighborhoods approximately neutral. Thus such configurations 
require the possibility of specific foldings of the PA at both local and global levels. 

The primary focus of the DH-type approach is the minimization of the extensive part 
of the energy by creating a homogeneous liquid-like structure, while non-extensive energies 
due to surface tension and electrostatics are relegated the role of determining the overall 
shape of the globule. By contrast, an RG-inspired approach |l(J to the problem assumes 
that the blobs form a self-similar structure (as depicted in Fig. |l|b) which attempts to take 
care of energies on every length-scale. In d space dimensions, the dimensionless interaction 
parameter at scale £ is u(£) = q\ij {Tr'jr 2 ). This expression represents the typical interaction 
energy of a random ^-monomer segment, assuming that the Coulomb interactions can not 
be screened. For a self-similar structure with oc t v ' , u(£) oc £ 1 ~^ d ~' 2 ' u , an d the interaction 
parameter grows or diminishes under rescaling unless v' = l/(d — 2). (The analogous 
argument leads to the exact value of Vh = 2/(d — 2) for uniform polyelectrolytes where no 
screening is possible ||20|| .) This result is valid only for 3 < d < 4: For d > 4 the electrostatic 
interactions are irrelevant and v = 1/2, while for d < 3 the polymer is stretched [y 1 = 1). 
Keeping interactions equally strong on all length-scales also generates a condensation energy 
propotional to Nq^/a. However, as this argument does not provide the prefactor, i.e. the 
actual value of the condensation energy, it is not possible to deduce whether the DH-type 
or RG-type ansatz produces the lower free energy state. 

Are the two aforementioned approaches mutually exclusive? In this work we present 
evidence that they actually represent two facets of the same problem: DH-theory attempts 
to minimize the condensation energy without paying attention to the surface. However, 
an object can be (locally) "compact" and still have an extremely extended shape which is 
controlled by the non-extensive part of the energy. The RG-inspired approach attempts to 
accommodate the latter energy. 

III. NUMERICAL SIMULATIONS 

Configurations of a polymer are completely specified by listing the position verctors {r^} 
(i = 1, . . . , N) of its monomers. The shape and spatial extent of the polymer are roughly 
characterized by the shape tensor, 
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1 N i N N 

= X/ r «M r ^ _ jy2 X/ r *M X/ T 3 V ' (3) 
j=l i=l j=l 

where the greek indices denote the Cartesian components of the vectors. Thermal averages 
of the ordered eigenvalues Ai > A2 > A3 of this tensor (sometimes referred to as moments of 
inertia) are used to describe the mean size and shape; their sum, i.e. the thermal average of 
trS, is the squared radius of gyration R 2 . Since we are dealing with sequences of quenched 
disorder, these quantities must also be averaged over different realizations. In three dimen- 
sions, uniform uncharged polymers in good solvents are swollen; their R g scaling as N u with 
v = 0.588 as in self-avoiding walks. Polymers in poor solvents are "compact", i.e. described 
by v = §. 

The Monte Carlo procedure used in this work is identical to that described in Ref. JT3 



Here we describe some important features, while the more technical details are described 
in Appendix [A]. The simulated chains are composed of iV monomers whose positions are 
discretized to a cubic lattice (d = 3) with lattice constant a. The connectivity of the 
polymer is maintained by restricting the maximal distance between neighbors to 4a. The 
excluded volume interaction is enforced by not allowing two monomers to come any closer 
than y/2a. Each quench is characterized by a set of charges = ±g - The electrostatic 
interactions between the charges, U = Y^aj) £4/(l r i — r j\), is included by assigning energy 
Uij(r) = qiqj/y/c + r 2 to each pair at a separation distance r, with c = 2a 2 , which 
"softens" the potential at short distances. 

The results of the simulations are parametrized by the chain length N, temperature T, 
and the overall excess charge Q = Each Q can be obtained by many realizations 

of randomness, and all results were averaged over 10 different quenches. However, rather 
than taking the same configurations through changing temperatures, 10 distinct quenches 
were used for each T and Q. The smoothness in variations of various quantities with tem- 
perature then provides added confidence in the thermal and quench averaging process. Not 
surprisingly, as explained in Appendix [A], the overall uncertainties are entirely due to quench 
averaging as the statistical errors of the thermal averages are smaller than the differences 
between quenches. 

Fig. depicts the temperature dependence of R 2 for 64-monomer chains. The number 
near each curve indicates the charge, Q/qo- At very high temperatures the electrostatic 
interactions are unimportant and the chains behave as self-avoiding walks, with R g cx N u 
and v = 0.588. The typical electrostatic energy of such configurations is estimated as 

« 2 - 1 2 N) 



(c,)K 5«\^; Ki ^' (4) 

where we have employed q[qj = (Q 2 — qQN)/N 2 for i j, and used R g as a measure of 
interparticle separation. Note that the interaction changes sign at 



Q c = q VN . (5) 

This is because the energy of strongly charged polymers is dominated by the repulsive 
interaction of excess charges. However, for weakly charged polymers, there is an attractive 
interaction between fluctuations in the charge distribution; the typical fluctuation of Q leads 
to the above result. 
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As temperature of the chain is lowered, the effects of interactions become apparent for 



Tq « (U) « (Q 2 - q*N)/aN" . (6) 

Chains with charge larger than Q c expand, while those with Q < Q c shrink with decreasing 
temperature. The strongest deviation occurs for the fully charged polymer, which for N = 64 
starts at Tq ps 320g 2 /a. For Q = 0, the deviation begins at the much lower temperature 
of Tq pa 5g 2 /a. Indeed, on the logarithmic scale of Fig. |2| the departure from the infinite 
temperature values of R 2 is most apparent for Q = 64g , starting at ln(T) ~ 6 (beyond the 
limits of the figure), compared to ln(T) ~ 1.7 for Q = 0. 

To see if the averages in Fig. |2] provide a good measure of the PA size at low temper- 
atures, we constructed histograms of the distribution of R 2 at T = 0.05g 2 /a for several 
values of Q. As thermal fluctuations are small, the histograms in Fig. |3] represent differ- 
ences between quenches. The distributions are fairly narrow, their widths not exceeding the 
distance between their averages. Thus a point in Fig. § provides a good measure of R 2 (Q), 
independent of further details of the sequence. The average of R 2 for unrestricted quenches 
is then obtained from R 2 g {Q) as, 

poo 

^(random) = J dQR 2 g (Q)P(Q) , (7) 

where P(Q) oc exp[— Q 2 / (2q 2 N)] is the probability density of an excess charge Q. In pre- 
vious work we found a very broad distribution for R 2 (random). Even in sampling a 
few quenches, there were several completely collapsed, and some strongly stretched configu- 
rations. Fig. [| shows the spatial conformations of several quenches examined in this study. 
The weakly charged configurations for Q/qo = or 4 are spherical globules, indistinguishable 
from each other. The chains are slightly expanded for Q/qo = 8, while for Q/qo = 16, a value 
of only twice Q c , they are strongly stretched. It can now be appreciated that the previously 
observed breadth of the distribution for R 2 (random) simply follows from the strong depen- 
dence of R 2 on Q rather than indicating a large scatter of R 2 amongst different quenches 
with the same Q. 

The averaged radii of PAs in Fig. || change monotonically with temperature. This sug- 
gests that compact and extended states are separated in the (Q,T) plane by a straight line 
starting from Q = Q c = qoVN at infinite T. This hypothesis was tested by looking at 
the Q and N dependence of the radius of gyration for chains of lengths iV = 16, 32, 64, 
128. To achieve good thermal averages, simulations were performed at T = 0.1g 2 /a and 
not at the lowest temperature in Fig. The dependence of R 2 on Q is depicted in Fig. ||. 
The vertical axis is scaled by N 2 ^ 3 to remove the A-dependence of the R 2 of the compact 
globules at Q = 0. The charges on the horizontal axis are scaled by Q C (N) for all polymer 
lengths. Although monotonic, R 2 exhibits strong variations with Q. The radius is barely 
increasing for small Q, but an extremely steep rise begins beyond a threshold charge. Due 
to the monotonic increase or decrease of PA sizes with temperature, this variation becomes 
even sharper at low temperatures. Fig. ^ strongly suggests that the transition from compact 
to stretched configurations at low temperatures still occurs for Q w Q c . 

When the distortions of a globular object are small, the changes in R 2 are not very 
sensitive to the changes in shape. The increase in the largest eigenvalue of the shape tensor 
Ai due to elongation of the object is partially compensated by a decrease of other eigenvalues. 
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Fig. H is analogous to Fig. |^ except that it depicts the charge dependence of the Ai. The 
strong elongation of the PA for Q > Q c is also apparent from this figure. 

IV. ANALOGIES TO CHARGED DROPS 

To explain the above results, we start with the empirical observation that PAs with 
vanishing excess charge Q compactify to spherical "globules" of spatial extent R ~ aN 1 / 3 
and surface area S ~ a 2 N 2 / 3 . It is thus natural to represent the total energy (or rather the 
quench-averaged free energy) of such globules as a sum of condensation and surface energies, 

E PA (Q = 0) = -e c N + 1 S . (8) 

The condensation energy is proportional to q 2 /a, while the surface tension is 7 = pq^/a 3 , 
where the dimensionless prefactor p « 0.1, is found to be rather small [^TJ. It should be 
emphasized that Eq.@ is not self-evident as it represents the average energy of a connected 
chain of N monomers with long range interactions, rather than N independent particles. 
While the existence of an extensive condensation energy is natural, non-extensive corrections 
may in principle be present without any relation to the surface. The presense of the surface 
term is deduced from the numerical observations that the object is approximately spherical. 
We may hope that this form of the energy persists as long as the deformations of the globule 
are not too large. 

If we now uniformly add a very small number (Q/qo <C A 1 / 3 ) of charges along the chain 
(e.g. by randomly replacing Q/2q of negative charges by positive ones), without modifying 
the spatial conformation of the PA, its total energy increases approximately by qo\Q\/a, 
representing the sum of changes in local interactions. For moderate charges (Q > qoN 1 ^ 3 ) 
the energy increase is dominated by the long range interactions and is of order Q 2 /R. If the 
PA is now allowed to relax, it will lower its energy by finding more favorable configurations, 
and thus, 

E PA (Q) < ~e c N + 1 S + Q 2 /R . (9) 

The above considerations are equally applicable to a charged drop, and we shall explore 
such analogies to treat weakly distorted PAs. In the following paragraphs, we initially 
review the results pertaining to the shape of a charged conducting drop. This analogy is 
most appropriate for an annealed version of the problem in which the charges are free to 
move along the polymer chain. We then go on to consider the shape of a charged insulating 
drop of immobile charges, which is a better representation of quenched PAs. 

For a conducting drop all the charge accumulates on the surface. The non-extensive 
contribution to the energy of a freely suspended spherical drop of radius R is 

E{Q) = \^ + ^R 2 1 . (10) 

However, the drop can change shape to minimize the sum of surface tension and Coulomb 
energies. The surface energy of an uncharged drop, E(0) = 47iR 2/ y, sets the overall energy 
scale of the problem while the dimensionless parameter, 



8 



a = 37(167^7) ee g 2 /g^ , (11) 

determines its shape. We shall refer to Qr as the Rayleigh charge of the drop. Note that the 
estimates of 7 and R following Eq. @ for the model PA lead to Qr ~ Q c . Let us initially 
consider only small deformations in shape of a single drop. Investigations of the shape of 
charged liquid drops go back at least to the last century p2| |. Some of the several variants of 



this problem are: (a) a freely suspended charged conducting drop Pffiffil ; (b) an uncharged 



conducting drop in an external electric field p8p9| ; (c) a dielectric drop embedded in a 
different dielectric liquid in an external electric field fl30"H34f . The last problem has received 
recent attention because it is mathematically identical to that of a drop of magnetic fluid in 



a magnetic field ||32l,34j-|37|l. The conducting drop is the limiting case of an insulating drop 



with infinite dielectric constant. Problems (a) and (b) are strongly related since similar 
shape instabilities are induced by both the external field, and by the internally generated 
field of a charged drop. In fact, in many experimental situations (see, e.g. Ref. J5H]), the 
drop is suspended from a tube. Creating instabilities by raising the potential of the tube 
is intermediate between the idealized situations described by (a) and (b). The stability 
conditions in experiments are usually discussed from the point of view p3| of increasing 
potential rather than charge. Therefore, many features of the problem which are of interest 
to our study are usually not addressed. In the Appendix [B| we summarize the relevant 
aspects from our perspective, while only quoting the main results in this Section. 

A closely related problem is that of an insulating uniformly charged drop. This problem 
has been considered in the framework of the charged drop model of atomic nuclei (see, e.g., 
Ref. [P9"| ). The non-extensive portion of the energy of a spherical shape is now 

E^(Q) = ^ + 4n 1 R 2 . (12) 

The similarity between Eqs. (|TT]) and (0) is evident. As explained in the Appendix [FJ this 
similarity persists even for the non-spherical shapes discussed in this Section: all results for 
conducting drops are transformed into results for insulating drops by replacing Q 2 /2 with 
3Q 2 /5. This analogy, however, does no persist to arbitrary deformations and differences 
between conducting and insulating models will become apparent in Section [V]. 

If the only allowed deformations of the charged drop are to ellipsoids of rotation (prolate 
spheroids), the spherical shape remains stable until a reaches 0.899. At this point the 
drop becomes strongly elongated with eccentricity e = 0.95, and continues to stretch with 
increasing a. For large a the long axis of the drop (and hence R g ) is proportional to a 1 / 3 . 
Due to the sharp increase in the aspect ratio of the spheroid with increasing a, the increase 
in the energy of the system slows down and becomes of order of (a In a) 1 / 3 (compared to 
order of a for the undistorted sphere). Fig. |^ depicts the resulting dependences of R 2 g and E 
on a. The behavior of R 2 g in Fig. [7| closely resembles the sudden expansion of polyamphilic 
gels in ref. 0. 

A more quantitative comparison between our results and the predictions of the charged 
drop model is possible: The transition in Fig. [5] appears at a' = Q 2 /(Qo^0 ~ !• The 
instability of a charged drop to ellipsoidal shape occurs for a ~ 0.9 in the conducting case, 
and a ~ 0.7 in the uniformly charged case. In Eq. |BT] we show that a = Q 2 / (12 V7), where 
V = a 3 N is the volume of the system, while 7 ~ 0.1q 2 /a 3 [^TJ. Thus, a ~ a', and the 
observed critical value of a' is surprisingly close to the predictions of the model. (Given the 
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numerous approximations of the model, such excellent agreement is probably fortuitous.) 
We conclude that, as long as the value of a is not too large, our MC results, the predictions 
of the spheroidal drop model, and experiments, are in good agreement. 

For comparison with Fig. [7], Fig. |8] depicts the energy per monomer as a function of 
scaled excess charge, obtained from MC simulations. At Q — 0, the curves for different As 
almost coincide. (Slight corrections of order A -1 / 3 are present, but invisible at this scale.) 
For small Q, the energy per monomer increases as Q 2 /{RN) ~ Q 2 /A 4 / 3 , as indicated by the 
straight dashed line in Fig. ||. For large charges, the energy increase slows down, indicating 
distorted PAs. Since distortions start for Q 2 /q%N « 1, in terms of variable used in Fig. |8| 
the departures of longer chains begin earlier. 

The results of this Section are not sensitive to the exact shape of the elongated drop. 
Any shape characterized by a long dimension R\\, and a short dimension R±, (such as 
major and minor semi-axes of an ellipsoid), reproduces the same answers qualitatively. The 
electrostatic energy is approximately ~ Q 2 /R\\, while surface energy grows as jR±R«. Their 
sum has to be minimized subject to the constraint of fixed volume, imposed by requiring 
V ~ R\R\\. For example, Gutin and Shakhnovich |5J consider the more general case of 
Q ~ N@. Minimizing the total energy for elongated shapes, they find Ru ~ A^ 4 ^ 1 " 3 , and 
R± ~ N 2 ^-/ 3 )/ 3 . Only a finite stretching is predicted for (3 = 1/2. (A directed version of this 
problem also exhibits a continuously varying exponent u((3) ^IJ.) Another recent study by 
Dobrynin and Rubinstein ||2|| relaxes the constant volume constraint and also reaches the 



conclusion that there is an onset of stretching for (3 = 1/2, although a completely stretched 
state is reached only for (3 = 2/3, when the Coulomb energy becomes extensive. 

However, as we shall show in the following sections, the ground state of a charged drop is 
not a simple elongated shape for large values of a. Conducting drops can shed away excess 
charge, while insulating drops disintegrate in a process similar to nuclear fission. Related 
pathways are available to PAs. 



V. BEYOND SINGLE DROPS 



Linear stability analysis indicates that a spherical shape is unstable to a variety of per- 
turbations |p6| |. Experiments show that a conducting drop disintegrates when the Rayleigh 



stability limit is exceeded. Nuclear fission demonstrates the corresponding instability of 
insulating drops. Appendix presents several mechanisms by which conducting and insu- 
lating drops can decrease their energy. For example, we show that a conducting charged drop 
can get rid of its entire electrostatic energy by emitting an infinite number of infinitesimal 
droplets. 

The mechanisms discussed in Appendix y rely on the breakup of the charged drop. 
Such routes are not available to the PA chain which must maintain its connectivity. Is 
the charged PA susceptible to similar instabilities despite its connectivity? Fig. |9] depicts 
Rg(T = O.lgo/ 'a) / R? g (T = oo) as the function of the reduced charge Q/N 1 ^ 2 for different 
values of N. The curves become steeper with increasing N and intersect at Q/N 1 ^ 2 « 1.4g - 
At the intersection point the radii scale as self-avoiding walks (with a prefactor slightly 
larger than the infinite temperature value. For Q > 1AQ C the PAs at low T are more 
stretched than self-avoiding walks, in disagreement with the finite elongation predicted for 
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ellipsoidal shapes. How does the PA go beyond the ellipsoidal limit while maintaining its 
connectivity? We can still exploit analogies to charged drops. 

The annealed PA is modeled by a deformable conducting "drop" of particles of size a, 
constrained to maintain a spherical topology. Although such a drop cannot expel charged 
particles, it can still reduce its energy by extruding charge in a finger of length L and diameter 
a. Balancing the Coulomb energy (Q 2 /L) of the finger with the increase of surface energy 
7aL, we find that the optimal finger length is L « Q/^ja. Fingers appear spontaneously 
only if their cost (roughly Q^ya) is less than the Coulomb energy of the uniformly charged 
sphere, Q 2 /R, i.e. for Q > R*/ja. The fingering instability occurs for a ~ a/R -C 1, i.e. 
far below the shape instability of a sphere. Thus the typical annealed PA has a protruding 
finger of length L oc Q oc N 13 for (3 > 1/3. But, as the weight of the finger is small, it 
does not effect the scaling of R 2 g for (3 < 5/9. Such PAs have large spanning sizes without 
appreciably greater R g . 

The insulating (uniformly charged) drop follows a different route. As is shown in Ap- 
pendix 0, it is stable for small a, but reduces its energy by splitting into several droplets 
of equal size for a > 0.293. We can again constrain the overall object to remain singly 
connected by linking the droplets via narrow tubes of total length L and diameter a. As 
long as La 2 <C R 3 , most of the charge remains in the spheres. The total electrostatic energy 
is proportional to Q 2 /L, while the surface energy cost grows as jaL. Equating the two gives 
L oc Q; not surprisingly, of the same order as the fingers in the conducting case. However, 
whereas the surface tension in the conducting case results in one big central drop, for the 
insulating case the droplets are separated as in a necklace. The radius of gyration is now of 
the same order as the span of the necklace, L. 

VI. RANDOMNESS IN THE NECKLACE MODEL 

The necklace model provides a good picture of a polymer with a short-range attractions 
between its monomers and a uniformly distributed excess charge: A polymer with Q 3> Qr 
is split into roughly a oc Q 2 /N beads connected by a string. Each bead is just below 
the Rayleigh threshold, and the string is stretched by their Coulomb repulsion to a length 
L oc Q. (Note that, as shown in Appendix [C], the optimal number of beads is proportional to 
(Q/Qr) 2 and not Q /Qr.) For Q <C qoN, only an infinitesimal fraction of monomers are part 
of the string, and the overall extensive part of the energy is unchanged. This picture should 
extend to any deterministic sequence, e.g. composed of alternating charges WJl, which has 
a compact state when uncharged. 

Is the necklace model also applicable in the presence of random charges? For our model 
PAs a m Q 2 / (q^N), with a prefactor almost identical to unity, and we shall use this relation 
as an exact definition of a. For a > 1 we may try to split a chain into a segments of 
approximately equal size. Each segment has average charge Q/a oc N/Q and incorporates 
N/a oc (N/Q) 2 monomers. Thus the fluctuations in charge of each segment are of the order 
of the average charge itself, and the picture of uniform, mutually repelling, beads is no longer 
applicable. It is not clear how we should model the shapes and distribution of the segments 
which have as of order one. 

Let us illustrate the difficulties caused by randomness for the case of an unrestricted PA. 
Since Q 2 = q$N, where the overline denotes an average over the ensemble of all quenches, 
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we have a = 1. As demonstrated in Appendix |C], the insulating drop is unstable to splitting 
already for a ~ 0.3, and thus a typical random PA is expected to form several globules 
connected by narrow tubes. Now consider splitting sequences of N monomers with total 
charge constrained to a particular a into two equal subchains of charges Qi and Q2- It is 
easy to show that each segment has a su bchain = (1 + ot)/2, while the mean product of the 
charges is Q1Q2 = <?oA(a: — l)/4. The subchains have, on average, values of a close to unity. 
Also, for a — 1, the average value of the product of charges vanishes. We thus have the 
paradoxical situation in which most spherical shapes are unstable, while there is on average 
no energetic gain in splitting the sphere into two parts. It is most likely that the ensemble 
of chains with a ~ 1 contains a broad distribution of sizes and shapes. 

Thus charge inhomogeneities drastically modify the necklace picture. The resulting PA 
is probably still composed of rather compact globules connected by a (not necessarily linear) 
network of tubes. The globules are selected preferentially from segments of the chain that 
are approximately neutral (or at least below the instability threshold), while the tubes are 
from subsequences with larger than average excess charge. It is amusing to inquire how a 
random sequence is best partitioned into large neutral segments. The resulting segments 
appear to have a broad distribution which will be addressed in future publication . 

In summary, we find that the behavior of PAs, and other charged polymers is controlled 
by the parameter a cx Q 2 /N. Chains with small values of a form compact spherical globules. 
The globules split for a > 0.3, resulting in a necklace of beads if the charge inhomogeneity 
is small. The span of the uniform necklace scales with the net charge Q. We don't have 
a consistent theoretical picture for the random PA beyond the instability threshold. The 
numerical results in Fig. |9] suggest that the size of such PAs grows faster than that of a self- 
avoiding walk, i.e. v > 0.6. The simulations so far are not inconsistent with v = 1 suggested 
by a scaling argument []10[ . However, as the simulations suffer from the usual shortcomings 
of small sizes, sampling, and equilibration, a definitive answer about the behavior of PAs is 
still lacking. 
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APPENDIX A: MONTE CARLO PROCEDURE 



The model used for Monte Carlo (MC) simulations is described in the Section [D J 



Dis- 

cretizing monomer locations simplifies checking for excluded volume interactions; allow- 
ing the bond length between the nearest neighbors to fluctuate without energetic cost (a 
"square-well" potential) provides sufficient flexibility to facilitate equilibration. The square- 
well potential has been used before in continuum simulations of tethered surfaces [ 43 1 ; on 
discrete lattices it is known as the fluctuating bond method [54 . 
procedure are as follows. 



The details of the MC 
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For each T and Q we start by selecting a (quenched) random sequence of charges ±go, 
whose sum (the total excess charge) is fixed to Q. This is accomplished by randomly selecting 
(N — Q/q )/2 positions on the chain for negative charges, and placing positive charges on 
the remainder. For each quench, we perform a thermal equilibration at temperature T, 
and then calculate the thermal averages of interest. The thermalization is repeated for ten 
different sequences and the results are averaged over the quenches. In an elementary MC 
step a monomer is picked at random and moved a single lattice unit. The move is accepted 
according to the usual Metropolis rule. Since each move requires recalculation of interaction 
energies, it involves O(N) operations. The MC time unit is defined as the period during 
which iV attempts are made. Thus, the CPU time per single MC time unit increases as N 2 . 

Obtaining good averages in random systems is a significant challenge. Errors appear 
due to both inadequate thermal equilibration and insufficient quench averaging. Our high- 
temperature chains (for T > hq^/a) are similar to uncharged polymers and their equilibration 
is limited by the slowly decaying "Rouse modes" . The slowest decay time is approximately 
the interval taken by a polymer to diffuse its own radius of gyration, estimated as follows: 
Since the acceptance rate of an elementary MC move is of order one throughout the simu- 
lation, the diffusion constant of a single monomer is also of order one (in units of squared 
lattice constant divided by the MC time unit). The diffusivity of the polymer center of 
mass is N times slower, resulting in a diffusion constant of D pa a 2 /N, and a relaxation 
time of t' = R 2 g N/a 2 . At high temperatures r' scales as N 1+2u , where v = 0.588. At low 
temperatures (T < OAq^/a) the polymer is almost compact, and a characteristic time can 
be obtained by considering phonons, plasma oscillations, or large-scale density fluctuations. 
Such time scales, in our MC time units, grow as r" ~ R 2 /a, 2 ~ N 2 ^ 3 . Unfortunately, there 
are probably much slower (and more important) time scales associated with crossing over 
large barriers to shape rearrangement which are thermally activated. We have no estimates 
for such times. 

We used r = iV 2 MC units as the basic equilibration time. Each equilibration lasted 
250t, but the first lOr configurations were dismissed in calculating thermal averages. As 
the number of operations per equilibration increases as N A , this is close to the maximal 
equilibration time which can be reasonably used in a simulation of this type. Several hours 
of CPU time (on Silicon Graphics R4000 workstation) were needed to equilibrate each quench 
at a given temperature for N = 64. Consequently, more than a day of CPU time is used 
to generate a single date point by averaging over 10 quenches. To collect all the data on 
N = 64 chains we needed about two months of CPU. For N = 128 we spent 10 days of CPU 
to obtain a single data point, and therefore only the Q dependence at a single temperature 
was investigated. 

We believe that the times used in equilibration produce satisfactory thermal averages. 
A direct check of the temporal correlation function of the radius of gyration for N = 64 
indeed indicates that the correlation time is approximately equal to r at high temperatures. 
This suffices to produce very good thermal averages. For example, a particular sequence of 
N = 64 monomers with Q = at T = 25 has average R 2 = 152a 2 , with standard deviation 
of approximately 60a 2 . For this polymer r' = 152 * 64 = 2.5r, and thus our simulation 
contains approximately 100 independent configurations. Therefore the average value of R 2 
is accurate to about ±6a 2 . However, the average R 2 for 10 distinct quenches are scattered 
over an interval of width 40a 2 . Thus the accuracy of our thermal averaging suffices to show 
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that different quenches have slightly different thermal averages of R 2 g . The accuracy of the 
final average is, therefore, limited by the number of quenches rather than by the thermal 
averaging. 

The correlation time obtained at low temperatures from temporal correlation functions 
of neutral PAs is shorter. This just reflects the reduction in size of the entire polymer. 
As mentioned earlier, correlations of R g are rather insensitive to shape changes and their 
much longer activated time scales. To obtain some, admittedly indirect, measure of the 
quality of equilibration for dense polymers, we compare our simulations with the quite 
extensively investigated restricted primitive model (RPM). The latter represents a solution 
of positive and negatively charged particles (±<7o)) interacting via a Coulomb force and a 
hard core repulsive potential of diameter a. (For a review of the subject see Ref. [(45H .) 
The thermodynamics of the model is conveniently presented in terms of a dimensionless 
density p* = na 3 where n is the actual number density, and a dimensionless temperature 
T* = k B Tajq\. At low temperatures the solution undergoes a phase separation between 
high and low density phases. As indicated by the dashed line in Fig. [1^, the critical point 
occurs at very low density (unlike regular fluids with short range interactions). Numerical 
investigations of this phase transition |^6|-^9[1 have encountered considerable difficulties: 
despite the low density, the behavior of the system becomes very erratic close to the critical 
point. 

Since our simulations also involve low temperatures and relatively high densities, it is 
interesting to find out where our system is located on the (T*, p*) -plane. Of course, as we 
are dealing with a single polymer rather than a dense solution, the comparison involves few 
somewhat arbitrary factors. Our lattice potentials approximately mimic interactions of the 
hard core particles of RPM. From this comparison we relate the MC temperature to T* by 
T* = \.2Tajq\. Secondly, we calculate the polymer density, assuming that the monomers 
uniformly occupy the volume of a homogeneous ellipsoid with identical eigenvalues of the 
shape tensor {Aj}. This leads to a reduced density, p* = 2.4a 3 /v / AiA 2 A 3 . For a neutral 
PA, p* is approximately independent of T for T > 5q^/a, and increases at lower T, leading 



to the trajectory indicated by the solid line in Fig. [10]. At densities close to the critical 
density of RPM, the PA temperature is almost an order of magnitude higher than the critical 
temperature. At lower temperatures the trajectory of PA simulations approaches the phase 
boundary on the "liquid" side, at densities twice higher than the critical density of RPM. 
Thus our polymers stay away from the problematic region where critical fluctuations may 
cause significant equilibration problems. Since our equilibration times exceed by several 
orders of magnitude those used in RPM model simulations, we believe that we have well 
equilibrated results. 

The acceptance rate of MC moves in our simulations is approximately 0.6 for all tem- 
peratures. It drops to 0.46 at T = 0.05<2o/a, and further lowering of temperature leads to 
a gradual "freezing". Repeated heating and cooling cycles performed on several samples 
indicates that the behavior is essentially reversible for T > 0.05gQ/a. We believe that the 
configurations obtained for T = 0.05<2o/a are very close to the actual ground state. Most 
of the low temperature investigations were actually performed at T = OAq^/a where the 
thermal averages are more reliable. For chain length N = 64 we performed an extensive 
study of the dependence of R g and other quantities on T and Q. We also investigated the 
iV-dependence of these quantities at T = OAq^/a for iV = 16, 32, 64, and 128. 
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APPENDIX B: SPHEROIDAL DISTORTIONS OF CHARGED DROPS 



In this Appendix we discuss the minimum energy shape of a single charged drop. We 
first examine the conducting drop, and then relate the results to insulating ones in the last 
paragraph. In terms of the dimensionless parameter a, Eq. (|TTD for the energy of a charged 
conducting spherical drop is 

E(Q) = E(0)(l + 2a). (Bl) 

For large Q the spherical shape is unstable, and the stable shape is determined by minimizing 
the sum of surface and electrostatic energies. In an equilibrium shape the surface is an 
equipotential, since otherwise energy can be reduced by redistributing the surface charge. 
The pressure difference between the inside and outside of the drop at any point on the 
surface is given by, 

Ap = 7 ( — + —) - 2ncr 2 , (B2) 

where r\ and r 2 are the principal radii of curvature, and a is the surface charge density. 
In equilibrium the pressure inside the drop must be constant. Thus, Eq. ([B2]), with a 
constant value of Ap at all points, determines the equilibrium shape. This is in fact a rather 
complicated integro-differential equation whose general solutions are not known. Note that 
for Q = Qfi, we find Ap = for a spherical drop. 

In 1882 Lord Rayleigh investigated the stability of a charged drop p2] and showed that for 
a = 1 the sphere becomes unstable to surface distortions described by the Legendre function 
P2(cos#). (The points of instability for higher harmonics are given by a n = (n + 2)/4.) The 
instability does not result in a small distortion as a exceeds unity but, rather, leads to 
a strongly elongated shape. For a > 1 no exact analytical treatment is available. Some 
progress is possible by assuming that the drop is an ellipsoid of revolution (i.e. a prolate 
spheroid) PBfl . Since both the surface area and the electrostatic energy of such shapes are 
known (see, e.g., Ref. 150]), the problem reduces to the minimization of 



E(Q) - ^ 




+ 2a K - '— In ((1 + e)/(l - e)) 



(B3) 



with respect to eccentricity e = y 1 — b 2 /a 2 , where a and b are the major and minor semi- 
axes. The first term in the brackets is the total surface area (incorporating the constraint 
of fixed volume), while the second term is the electrostatic energy (Q enters via a). It has 
been shown [E^] that a prolate spheroid is not an exact equilibrium shape: for aspect ratio 
a : b = 2 : 1 the pressure difference Ap in Eq. (|B2| ) varies by about 2% for different points 
on the surface. Nevertheless, we may assume that the forces pulling the ellipsoid out of the 
shape are small as long as a/b is not excessively large. Numerical solutions of a drop in an 
external field confirm that the ellipsoidal approximation is reasonably good for shapes that 
are not too elongated p9||34| , |33| . 
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Ailam and Gallily |23] noted that Eq. ( P3|) has a local minimum for e ^ even for a 



smaller than unity, i.e. below the Rayleigh stability limit. However, they did not determine 
this range accurately. Fig. [11] depicts E(Q) as a function of e for several values of a in the 
relevant parameter range. A new local minimum first appears for a = 0.887, and becomes 
the global minimum for a > 0.899. At the latter a a spherical drop should discontinuously 
"jump" to a strongly elongated shape with e = 0.95. However, the spherical shape (e = 0) 
remains a local minimum until a — 1. In an ideal experiment in which the charge is gradually 
increased the drop stays in a metastable spherical shape. At a = 1, the sphere becomes 
unstable and stretches to e = 0.98. For a 3> 1, the eccentricity approaches unity. In this 
limit, the asymptotic forms of the surface and electrostatic energy are (1 — e) _1//6 and — (1 — 
e) 1 / 3 ln(l — e) respectively. The minimum energy is achieved at e 2 ~ 1 — (n/8a) 2 / In (ir/8a) 2 , 
while the energy increases as 

E(Q) « 2.03£(0)(odna) 1/3 . (B4) 

Note that the asymptotic increase of E(Q) (oc Q 2j/3 ) is much smaller than that of the 
undeformed sphere (oc Q 2 ). The exact dependence of E(Q) on a is plotted in Fig. [7]. In the 
same figure we also present the squared radius of gyration of the elliptical drop, 



R 2 9 (Q) _a 2 + 2b 2 _ (1 - e 2 )-! + 2(1 - e 2 ) 



R 2 (0) 3R 2 



(B5) 



Since for large Q, 1 — e 2 is proportional to a~ 2 up to logarithmic terms, R g ~ Rot 2 ! 3 ~ 
R(Q/Qr) 4 / 3 . While the ellipsoidal shape is not an exact solution to the charged drop 
problem, Eq. ( pi|) provides an upper bound on the total energy. 

We next consider a drop in which the charge is uniformly distributed over the volume. 
This system has been considered in the context of the liquid drop model of atomic nuclei 



see, e.g., Ref. p9|). The energy of the uniformly charged sphere, Eq. ([12]), is 



E®(Q) = E(0)(¥- a +l) , (B6) 
5 

As in the case of the conducting drop, the uniformly charged drop becomes locally unstable 
|5lH53l to infinitesimal distortions for a = 5/6. Assuming that the drop distorts into a 
prolate spheroid, its energy can be written down explicitly ||53|| . The resulting energy is 



identical to Eq. (|B3| ), except that the factor of 2a in the second term is replaced by 12a/5. 
The same factor relates Eqs. ( |B1| ) and (|B6|). Thus all the previous results for conducting 
drops are also applicable to insulating drop after multiplying as by 6/5. We should note 
that this simple substitution does not hold for drops of arbitrary shape. 

APPENDIX C: SPLITTING A CHARGED DROP 

An elongated ellipse is not a local equilibrium shape of a drop, since perturbations 
P n (cos8) with n > 2 become unstable pB| . There are other theoretical indications that 



no elongated shape is the global ground state. For example, as the eccentricity increases the 
electric field at a tip (~ Q/b 2 ~ Q(l — e 2 ) -1 / 3 ) becomes strong enough to support conical 
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tips [f23|,[33|]. Experimentally, it is observed that, for a > 1, a conducting charged drop disin- 
tegrates into smaller ones. Assuming that the experiments can be described by the charged 
drop model [p5| , p^ ,p5[] , disintegration begins for a equal to or slightly above one. This is not 
surprising, since already at a = 1 the equilibrium shape is strongly deformed. The instabil- 
ities lead to ejection of smaller droplets whose size distribution is believed to be controlled 
by hydrodynamic effects. Apparently, after a significant elongation, many roads towards 
decreasing energy open up and the choice is made by dynamical effects. In this Appendix 
we shall investigate the splittings of a drop into several spherical droplets removed to infinite 
separations. 

First consider the splitting of a single conducting drop of radius R into two (secondary) 
droplets with radii R\ and R2. Using Eq. (|B1|) for the energy of each droplet, the total 
energy of the infinitely separated pair is 

E 2 (Q) = ^ + ^ + ^i{R\ + Rl) ■ (Cl) 

The charges satisfy q\ + q2 = Q, while the radii are constrained by R\ + R\ = R 3 to preserve 
the total volume. Minimizing E2(Q) with respect to q%, while treating q2 as a dependent 
variable, gives 

B ' Q (C2) 



Substituting this result into Eq. ( |CI| ) we find 

£ 2 (Q) = £(0)(^ + £r/) , (C3) 

where the reduced radii rj = Rj/R satisfy the fixed volume constraint J2j r j — 1- (Note 
that E(0) denotes the surface energy of the original drop.) The stationary value of E 2 (Q) 
is found by solving dE 2 {Q) /dr\ = 0, and treating r 2 as a dependent variable. The resulting 
equation has several possible solutions, including 

rl = 4 = \. (C4) 

However, the symmetrical solution is a minimum only for a > 1. Thus strongly charged 
drops would prefer splitting into two equal droplets. A second solution, 



1.2 




(C5) 



exists only for a < 1, where it is an energy minimum, less than that of a single sphere. Thus 
a weakly charged drop can always reduce its energy by splitting into two unequal parts. 

Further breakup of the drop is possible, and we next consider splittings into n droplets. 
The solution to this problem is analogous to the previous case. The distribution of charge 
among the droplets is still given by Eq. (|C2|) and E n (Q) has exactly the same form as Eq. 
( p3| ) with j summed from 1 to n. The search for extrema of E n (Q) leads to several solutions, 
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classified by two sets of droplets. One set contains m small droplets of reduced radius a, 
while the remaining m — n have larger reduced radii b (b > a), such that 



ma 3 + (n 



m)b 3 - 
ab 



b (ma + 



n 



m)by 



a 



(06) 



The solution for a particular m exists only in a finite range of as, and finding which solution 
represents the global energy minimum is quite cumbersome. One limit, however, is easily 
examined: Consider splitting the drop into one large droplet and m = n — 1 smaller drops. 
In the limit of n — > oo, a? pa a/n 2 , and b pa 1; the total volume, area, and electrostatic 
energy of the charged droplets, vanishes, while they, nevertheless, carry all the charge of the 
system. Thus the energy of the system is reduced to E(0), i.e. the energy of the uncharged 
original drop! 

We next consider the insulating uniformly charged drop. If such a drop is split in two, 
the charge of each droplet will be proportional to its volume, i.e. = QRf/R 3 . Using 
the energy of a single sphere in Eq. (|B6|), the total energy of a pair of infinitely separated 
spherical droplets is obtained as 



E®(Q) 



E{0) 



3 3 



(C7) 



The charges {qj} and the reduced radii {rj} satisfy the same restrictions as in the conducting 



drop. For a < 1/6 the only extremum is a maximum at rf 



1/2. For a > 1/6, this 



point is a local minimum. However, only for a > 0.293 is the resulting energy lower than 
that of the original drop. A similar scenario is found for splitting the drop into n secondary 
droplets. Larger values of a are needed to stabilize solutions with higher n, and there is 
an optimal number of droplets (all of the same radius) for each a. Since the energy of an 
n-drop system is, 



n 



E(0) 



3/Q 
5 V n 
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n \ 
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+ 47T7 
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the optimal n (for large n) is found from dE^ jdn = as 

24 



n 



-a 



The total energy of the optimal configuration grows as 

^pLai(Q) = 2.53£(0)a 1 / 3 . 



(C9) 



(CIO) 



Thus the uniformly charged drop can not lower its energy as drastically as its conducting 
counterpart. Since we have not exhaustively searched for other configurations, the above 
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result should be regarded as an upper bound to the ground state energy. Note that this 
bound has the same scaling (up to logarithmic corrections) on a as that of the highly 
elongated spheroid. 
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FIGURES 



FIG. 1. Qualitative views of the spatial arrangement of 'blobs' in a PA with quenched ran- 
domness. Electrostatic interactions within each blob are smaller than /cbT. Lighter and darker 
shades of the spheres denote predominantly positively or negatively charged blobs. The DH view 
assumes that the blobs can rearrange in a pattern (a) where interactions are screened on long 
distances. According to an RG-inspired picture the blobs form a self-similar pattern (b) with the 
same interaction energy on all length scales. 

FIG. 2. R 2 (in units of a 2 ) as a function of T (in units of q^/a) for several values of the excess 
charge Q for a 64-monomer chain. Each point is an average over the 10 independent quenches 
used at each temperature. The numbers near each curve indicate Q/qo- 

FIG. 3. Histograms of the distribution of the (ten) values of R 2 (measured for N = 64 at 
T = 0.05gQ/a) for several charges Q/qo, indicated near the histograms. 

FIG. 4. Spatial conformations of 64-monomer PAs at T = 0. 05^0 /a, for values of Q/qo equal 
to (a) 0, (b) 4, (c) 8, and (d) 16. Dark and bright shades indicate opposite charges. The diameter 
of each sphere is about 0.4 of the actual excluded-volume range. 

FIG. 5. Scaled R 2 as a function of Q/qo for chain lengths N = 16 (open triangles), 32 (full 
triangles), 64 (open circles), and 128 (full circles). 

FIG. 6. Scaled largest eigenvalue, Ai, of the shape tensor as a function of charge Q/qo for 
several chain lengths. The symbols are the same as in Fig. || 

FIG. 7. Radius of gyration of the minimal energy spheroid in units of the radius of gyration of 
the undistorted sphere (left), and the energy of the spheroid scaled to that of the uncharged sphere 
(right), as a function of a. 

FIG. 8. Energy per monomer in units of q\/a versus scaled excess charge for several chain 
lengths at T = O-Obq^/a. The symbols are the same as in Fig. [|. The dashed line is the energy an 
undistorted PA. 

FIG. 9. Ratio between squared radii of gyration at T = OAq^/a and T = 00 as a function of 
scaled excess charge. The symbols are the same as in Fig. ||. 

FIG. 10. High and low density fluid ("liquid" and "gas") transition line (dotted) in the plane, 
(T*,/}*), where T* and are the reduced temperature and density (see text) of the Restricted 
Primitive Model. The solid line indicates the trajectory of the neutral PAs used in our MC 
simulations. 
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FIG. 11. Energy differences between drops of prolate spheroidal, and spherical, shape (normal- 
ized to the energy of spherical drop) as a function of squared eccentricity e 2 . The graphs correspond 
(from top to bottom) to a = 0.88, 0.89, 0.895, 0.898, 0.90, 0.91, and 0.92. 
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